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| Abstract 

> ■ 

qq ■ We develop a mathematically precise framework for the Casimir effect. Our 

working hypothesis, verified in the case of parallel plates, is that only the regulariza- 
tion-independent Ramanujan sum of a given asymptotic series contributes to the 
Casimir pressure. As an illustration, we treat two cases: parallel plates, identifying 
a previous cutoff-free version (by G. Scharf and W.W.) as a special case, and the 
sphere. We finally discuss the open problem of the Casimir force for the cube. We 
propose an Ansatz for the exterior force and argue why it may provide the exact 
(i^ • solution, as well as an explanation of the repulsive sign of the force. 

<D 

43 ■ PACS number: 12.20.Ds 

> : 

1 A General Framework for the Casimir Effect 

a. 

Significant progress on the Casimir effect from the experimental point of view oc- 
curred in recent times [|]]. In spite of that, several theoretical problems remain, such 
as a real explanation of the sign of the force in the case of compact regions. The 
situation is worse with regard to a mathematically precise framework for the effect, 
due to the cutoff (or regularization) dependence of the energy, a fact emphasized 
by C. R. Hagen in M and somewhat less emphatically by P. Candelas in ||. The 
physical reason why divergences occur is well understood Q and is that the bound- 
aries are treated by quantizing the radiation field with mode functions || which 
are adapted to the type of (classical) boundary conditions (b.c), e.g., Dirichlet or 
Neumann. However, real boundaries consist of electrons and ions and such b.c. are 
not justified except if the particles act collectively in an essentially classical manner 
||, which is a priori not the case 0, and our ignorance in dealing with this fact is 
signalled by divergences. 
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Divergences are, of course, well-known in field theory, but they arise here in a 
different way, as explained above. Mathematical physicists, and several theoretical 
physicists, agree that a mathematically precise framework to cope with these di- 
vergences would be conceptually useful. Such frameworks exist in field theory (see 
H and references given there). A cutoff- free or "finite" theory of the Casimir 
effect (in the spirit of [0]) was attempted by G. Scharf and one of us (W.W.) in ||. 
It requires, however, the use of periodic b.c, which are unphysical in the case of the 
electromagnetic field. 

In this paper we present a thorough derivation of the results first announced in 
||. We reconsider the problem introducing ab initio an ultraviolet cutoff (1/A). The 
Casimir energy (CE) E vac (A) would diverge if the limit A — > was taken, but we do 
not need to do so, because the Casimir pressure depends only on the A-independent 
term in the asymptotic expansion, which is the (RI) cutoff independent term of the 
Ramanujam sum of a divergent series. This idea is due to K. Dietz |TIJ. In section 
we show how the result of || for the parallel plates is recovered as a special case. 
Some of the ideas of || are also used and summarized below, for convenience. 

Following |lTJ, consider an electromagnetic field at T = enclosed in cavities of 
identical shape, but made of different materials, the latter providing natural cutoffs 
for the high-frequency spectrum of zero point modes. The vacuum energy is thus 
given by 

E yaa = ~Y t w«C«W, (!) 

a 

with C a (A) material dependent cutoff functions dependent on a variable A with 
dimensions of length, which we normalize by 

C a {A)\ A=0 = l. (2) 

Since E vac has dimension (length) -1 in natural units, it may be written as an (asymp- 
totic) series 

E vac = a L 3 A~ 4 + aiL 2 A~ 3 + a 2 LA~ 2 + a 3 A _1 + a 4 L~ x + a 5 L~ 2 A + . . . , (3) 

where L is a length characterizing the spatial extension of the cavity. Dietz conjec- 
tured |IJ that by a theorem of Ramanujan the A-independent term a^L~ l in ([3D is 
independent of the regularization (i.e., of the set {C a (A)} in (|l|)) provided (||]) holds. 
We shall return to this conjecture later. 

In this paper we consider as in || the prototypical example of massless scalar 
field confined in a compact region K (a compact manifold with boundary) - the 
modifications introduced by considering the full electromagnetic field will be men- 
tioned later. We show that the uj a in ([IJ) should be identified with the eigenvalues 
of the square root of the Laplace-Beltrami operator. This is not unexpected, be- 
cause the relativistic energy is \k\ = (fc 2 ) 1//2 , but it has important consequences for 
expansion ([]). Consider || the field A(x) quantized in infinite space 

1 f d 3 k r ^ -> -\ 

A(x) = r / -= a{k)e- ik x + a + {k)e ik - x ; (4) 



(2tt)5 J V2u; 
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[A-(x),A+(y)] = 1 D< i + \x-y) 



(27T) 



d 3 k 
2\k\ 



—ik-x 



Att 2 (xq — iO) 2 — x 1 



(5) 
(6) 



Time evolution is generated by the Hamiltonian H = J d 3 xH(x), whose density 
can be written in the form 



H(x) 



dA 
dxo 



-A 



d 2 A 
dxl 



(7) 



Normal ordering is defined in momentum space. In order to go over to a geometry 
with boundaries, we should formulate it in x-space by the point-splitting technique: 



OA 
dx 



dA(x) dA(y) 
lim : 1 ' yy > 

y^x dx dy 



(8) 



lim 

y 



i f dA(x) dA(y) 
x\ dx dy 



1 d 2 



Di + \x-y) 



i dxl 



Finally, 



H(x) 



lim 

y->x 



1 dA(x) dA(y) 

2 dx dy 



(9) 



Taking into account that real boundaries consist of electrons and ions and the field 
which interacts with them is quantized in infinite space, we consider (|9D to be the 
Hamiltonian density describing the field both free and with boundaries. In the latter 
case, however, the first two terms in @ must be defined in the adequate Fock space, 
i.e., the concrete representation of the field operator is dictated by the geometry. 
Consider a compact region K and Dirichlet b.c. A(x) = for x G dK. Then A(x) 
may be expanded as follows: 



A(x) = V [a n u n (x)e-^ X0 + a+u n (x)e^ x °] , 



(10) 



where u n are normalized real eigenfunctions of the Laplacian in K, satisfying Dirich- 
let or Neumann b.c. (discrete spectrum): 



Au n (x) = ulu n (x) . 



(11) 



The concrete Fock representation is now specified by regarding a+, a n as emission 
and absorption operators ([a n ,a+] = 5 nm ) and defining the vacuum by 

a n Vl = Vn. (12) 
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We thus find in this Fock representation: 



1 / d V 1 d 2 



+ \]im-^{D^(x-y)- D^(x,y)} 

I V~+X OXn I J 



(13a) 



where 



d k\x - y , x, y) = i 7^ 



u n (x u„\y e 



-iu> n (xo—yo) 



(13b) 



and the semicolons in (|13a|) denote normal ordering with respect to the new emission 
and absorption operators a+ and a n . Notice that is the solution of the wave 



equation ODq = with initial conditions 



D { + \+0,x) = -^^r—: 
4tt z x z + tu 

<9 ^ +) )(0,f) = ±5(x), 



(14a) 
(14b) 



Dq\+Q,x) is the Green's function of the square root of —A in infinite space ||. 
Similarly 



yo=%o 



(doD^) (x,y) 

D^\+0,x,y) 



^{x-y) ; 



-(-A K ) *{x-y) 



(15a) 
(15b) 



where Ak denotes the Laplacian on K, with Dirichlet or Neumann b.c. 

We now consider two types of cutoff functions, one of them general, satisfying 
(El), the other special, of type 



C a (A) = C{Aw a ) 



16a) 



satisfying 



C(0) = 1 . 



(16b) 



We shall also interested in a particular case of ( |16a|) , namely 

C(k) = e~ Ak k > . 



(17) 
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In terms of the special choice (p]), we may, by (^), (P^|)-( P^d]) , compute a regu- 
larized vacuum energy density H vac (x, A) in the following way: 



n J 



(18) 



As an aside, notice that (|TTD corresponds to ascribe a small imaginary part — iA to 
— 2/o = T , and thus represents a "natural" choice, akin to the principal value in 
distribution theory [12|]. For this special case (|Tj 



1 d 

where P, P satisfy the "heat equation" 

^ K -(-A s )^P(x,y;A) = 0, (19b) 
with the b.c. 

P{x,y;A) = if x or y E OK , (19c) 

in the case of Dirichlet b.c. 

There exist methods to compute the asymptotic expansion (in A) of the quantity 
in brackets in ( |19a| ) [13], which solve the problem in principle, but the actual form 
(H), with the given coefficients, depends on the details of the discrete (eigenvalue) 
spectrum of (—A) 2. 

Let now L be a linear dimension of the compact region K = Kl and M a linear 
dimension of a region Km of which Kl is a subset. Typically, if Ki is a cube of side 
L, Km is a cube of side M > L concentric with Kl, and similarly for a sphere or 
other manifolds. It is correct to impose the same b.c. (e.g. Dirichlet or Neumann) 
on Km in order to define the outer Casimir problem [|14], [11]]. In fact, previous work 
on the sphere using the Sommerfeld radiation condition was not correct, although 
the results were right, because it did not lead to real eigenvalues Define 

E vac (L, A, M) = EZ er (L, A) + E^(L, A, M) , (20) 



where 



and 



E^(L,A)= [ d 3 xH(x,A), (21) 
Jk l 

£° a u c tor (L,A,M) = / d 3 xH(x,A). (22) 

Jk m \Kt 
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As previously remarked, if Dirichlet b.c. are imposed on K^, H (respec. H) is the 
density fll8"D with the {u n } defined by Dirichlet b.c. imposed on K L (respec. K L and 
Km)- Analogous definitions hold for other b.c. (e.g. Neumann or mixed). If ([J), 
(0) is adopted, the second sum in (|20|) refers, then, to the modes u n corresponding 
to the solution of ([□]) in Km\Kl, with the above-mentioned b.c. Suppose that 
both E l ™ er (L, A) and E°£ er (L, A, M) have asymptotic series (§, and let E^ cr (L)(= 
a^ ner /L) and E°^: er (L,M) be the corresponding A-independent terms. Then the 
Casimir pressure on the boundary surface pc(L) (a measurable quantity) is defined 
by the thermodynamic formulae (zero absolute temperature): 

Pc(L)=j/S ncT (L)-p c ^(L), (23a) 



where the relative minus sign takes into account that p£ uter refers to a normal vec- 
tor pointing inwards towards Kl, while p™ ner refers to a normal vector pointing 
outwards, and 



Q pinner ( f \ 

PT C \L) = - ™ ; (23b) 

OVinncH-kJ 

jg-W—Jim af-'g'^ . (23c) 



and an important feature of the thermodynamic limit |16| is that the derivative in 
( |23c| ) is taken with M fixed, only L varies. 

It is essential that the CE be independent of the cutoff function C in ([!]) or ( |16a| ) 
provided it satisfies (Q) or ( |16b|) . As remarked in |10|| , a necessary condition for this 
regularization independence (RI) to hold is that (|3|) contain no logarithmic terms, 
because, otherwise, the 'A-independent term" is obviously ill-defined. For the cube 
there are no such terms in (|^), but such is not the case for the sphere; however such 
terms may be omitted in the case of the sphere because they cancel in the expression 
E™ ei (L, A) + E°£ eT (L,A,M) which have an asymptotic series (g) as M -> oo, so 
that for the sphere of radius a we have: 

q sphere sphere 

fc(a) = - w *V = °T^F ' (24) 

A full proof of RI is given in section 2 for parallel plates. We shall leave a more 
detailed discussion of higher dimensional cases JH]] to a further publication, but we 
wish to make a few important remarks: 

a) the A-independent term in (§) should coincide with the A-independent term 
of the Ramanujan sum of a divergent series of positive terms, such as ([I]), with 
C Q (A) = 1 see |l8j (p. 318 ff.) and section 1. According to this concept, for 
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instance 



1 + 1 + 1-1 + 



2 (»,o); 



1+2+3+ 
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(ft,0) 



taking the origin as reference point (see |Tj|, 13.10.11). This is proved in 



section |2] for parallel plates, and is the basis of RI; 

b) the present definition of the CE is mathematically rigorous. In particular, the 
limit A — > is never taken. In fact, (0) shows that, in general, it does not 
exist (an exception is the Casimir effect for parallel plates with periodic b.c, 
see || and section 0). The reason for this is that we do not know how to 
treat the surface properly in microscopic terms, a formidable problem (see the 
conclusion) ; 

c) RI justifies the definition of the A-independent term in (fj) as the CE physically: 
it reflects the field theoretic structure of the vacuum state which is independent 



of the cavity materials [10 . It is also expected to be the only term in (m which 



contributes to the pressure: this was proved in |K| for parallel plates. 

Section |2| entails a complete proof of RI for parallel plates (a "limit" of a compact 
region), as well as the explanation of the "theory without cutoffs" for the case of 
periodic boundary conditions in ref . || . The sphere is also treated a an illustration, 
in section ||]. In contrast to the parallel plates, the force for the sphere is repulsive, as 
known since the pioneering work of Boyer [TO . The nature of our derivation, which 



is a adaptation of the ideas and results of |L5| and |14j to our framework, does not 
convey an intuitive "explanation" of the sign of the force. This is a difficult problem 
because the CE is a sum of fluctuations of the electric and magnetic fiels in the 
vacuum state. A basic issue is: if the flat parallel plate geometry is changed to a 
compact manifold with boundary, how does the sign of the force change and why? 
This question is most clearly analysed in the case of the cube, which is the simplest 
deformation of the parallel plates geometry. In section |] we consider the interior 
problem for the cube, using the method of the Poisson summation formula used in 
||. This method had already been used for the same purpose in [19]. Since this 
reference is not readily available, we include our (independent) derivation which 
generalizes [|1^] in the sense that we obtain the full asymptotic formula, and fits 
nicely into the present framework. It should be also remarked that it coincides with 
the numerical result of [Z0||. The inner problem leads, however, to an attractive force, 



while the result for the sphere leads us to expect a repulsive force. Therefore, the 
repulsive nature must be due entirely to the exterior pressure. 

In section [5| we introduce an Ansatz to solve the exterior problem for the cube, 
which leads to a repulsive force. We also have applied the Ansatz to the only known 
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soluble case with flat geometry, i.e., the case of parallel plates (Appendix A). From 
the analysis of this soluble case we identify the physical reason why our Ansatz does 
not modify the pressure: the Ansatz introduces some extra stresses, but these are 
parallel to the plates (faces of the cube) so that the pressure on the plates (faces) 
is insensitive to these additional stresses, then providing the correct result for the 
plates and (we believe) for the cube. If the latter conjecture is true, an "explanation" 
of the sign of the force also follows. This is left to the conclusion and open problems 
in section || 

We have used a very general class of cutoffs in momentum space for parallel 
plates. For the cube the proof is essentially the same as for parallel plates, but there 
are some subtleties in the case of the sphere which have not yet been fully worked 
out. Nevertheless, it is an open problem whether only the cutoff-independent part 
of the CE is relevant to the pressure, except in the explicit case of parallel plates 
r0 |. We shall admit this as a working hypothesis throughout. 



2 Parallel Plates 

We consider the problem of parallel plates, with distance d along the z-axis; take 
the positions of the plates at z = and z = d, and adopt the form (]16a|) in (fHf), 
Q2"0D with Dirichlet b.c. (Neumann b.c. yield the same results). The inner Casimir 
problem corresponds to the region Kl = Kd = {x G 1R 2 x [0, d]}, and the outer 
one to the region K R \K L = {x G R 2 x [d, d + R)} U {x G R 2 x [-R,0]}. The 
eigenfunctions associated to the inner problem are 

< ner (fc x , k y ) = —\R sin (^-z) e ^ x+k ^ n = 1,2, 3- • • , (25) 



2tt V d V d 
corresponding to the eigenvalues of (—A) 2 given by 



<EX = v(t) 2 + A; ' + ^' (26) 



in 



Tip . The outer eigenfunctions are 



u^\k x , ky) = ± J I sin (^(z - d)) e^ + ^) ; 



(27) 

uT ei ' 2 (k x ,k y ) = ^|sin (J£z) e^ +k ^ 



with eigenvalues 



<£\ = \I(Z) +kl + k 2 . (28) 



V RJ 



S 



We first adopt the choice QT7|). Introducing polar coordinates in the x-y plane, we 
calculate the first (inner) sum in ( f20|) (we do not integrate along (x, y) G M 2 , which 
would yield +00). The proper way to do this is to limit the (x-y)-jAane integration 
to a finite region with area A, and then take the limit for £ = % (this procedure 
yields the same results presented here and we omit it for brevity): 



2(2vr)< 



2d / dkk 6 e 



3 „-Ak 



(29) 



+ 2vr V / dfc fc 



Performing the change of variable fc^ = y i^j) + k 2 in the second integral in 
the r.h.s. of (|29|) we obtain 



2(2tt) 2 
d 



-2d I dk k 3 e- Ak 
d 1 



OO /»00 

+ 2ttJ2 dk' n 

I TiTT 

n=l J ~ 



e -Ak' n 



-6A- 4 + 



<9A 2 



1 4ze 
1 d 

A 2 e4r-1 



(30) 



we now use the expansion ([|TH|, p. 320) in (|3"0| ) 

OO 

_^ = i_U + Y(-i) k - l B k 

e* - 1 2 ^ V ; 



fe=l 



(2*)! 



(31) 



obtaining (i? 2 



and thus 



30^' 



^c er (A,rf) 



1 



1 7T 



rmner 
^Casimir 



4ttA 3 2 720c/ 3 

1 7T 2 



+ 0(A) 



(32) 



(33) 



2 720c/ 3 ' 

Two remarks are in order. The surface term — 4^3 in (|32"D are absent for periodic 
b.c, because the latter allow the term n = in ( |2"9"D which exactly cancels it. This 
explains the result of ||. The external CE is zero due to (p7|), (^8|) because, for the 
outer problem, d in ( |3"2"D is replaced by R, and thus in the limit R — > 00 



pouter ri 

'-'Casimir u " 



Finally, 



^Casimir 



1 7T 2 

2 720c/ 3 ' 



(34) 
(35) 
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The above energy is one half of the result for the electromagnetic field, due to the 
summation over the two polarization states in the latter. Notice also that, in natural 
units, S is of order (length)" 3 . 

An amusing aspect of the present derivation is that it seems to depend on the 
choice (|17D , i.e., of an exponential cutoff in ([2^) and (|30"D , which, due to (|3~ID , leads 
to (|32|). Consider now a general cutoff function (|16a| ). Omitting the volume term in 
we may write 



1 

Cac er (A,rf)= lim — (36) 

where 



n~»oo g7T 

m=l 







/m7T\ z „ / . / /rmrY- 



g{m) = / du \/u+ C AWn+ V/y 



( m7r \ 
d J 

It is of interest to compute 



/oo 
C(AVu) . (37) 



nw \ /m7r\ 2 ^ / , m7r\ 

<?«(m) = -2(— ) ^(a^-J ; (38a) 



,m W = -4-(_) C (a_)-2(-) - C« (A-) ; (38b) 



By [[18] (p. 326), under the following conditions ( |44"D and (flop on C: 



where 



and 



n 2d r°° 1 

9(m) - — / rfg g 3 C(A g ) + -0(0) — > S fc , (39) 

7T Jq Z n— >oo 

m=l 



- (2]feT2)I / VW% (2fc+2) (*M , (40) 



o 



^fc( a; ) — 4>k{ x ) mod 1 (i.e., equal to 4>k(%) for (41) 

< x < 1 with period 1) , 
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and 4>k are defined by 



! n! 



and 



^(o) = E(- i r 1 ^ (2r '" ) (°)- («) 

We changed the notation of |Jl8| : the Ck on pg. 326 corresponds to our £ fe . Notice 



that the second term in (|39|) corresponds to the subtraction of the vacuum term, 
which appears in a natural way as a necessary requirement in a purely mathemat- 
ical context! The term |g(0) contributes only to the A-dependent terms in the 
asymptotic series. 

Theorem Let the special cutoff junction oj 'type ( 1 6'a\ ) satisfy, besides ( \16't\ ), the 



conditions: C is infinitely differentiable and its derivatives (C^ = C) satisfy 

/oo 
C {k \x)dx < oo ; (44) 

C {k \x) — ► . (45) 

x— >oo 

Then, for Dirichlet (or Neumann) b.c. the A-independent term in (Qj is the cutoff- 
independent part of the Ramanujan sum of the divergent series ffi) with C a (A) = I, 
and is therefore RI, i.e., independent of C. 

Remark (k > 1) is referred to as the (9?, 0) sum of the (divergent) series 
Y2m=i 9i m )i where 3? refers to Ramanujan and to the reference point (the origin in 
our case). Usually (see, e.g., f2~l"fl , p. 138), the result is presented informally without 
the important last term in (flU)), and assuming that C satisfies C (fc) (0) = for all 
k > 1, besides ( |16b| ), which is not satisfied by the special choice ( |TTD (see, however, 



22| for a nice approach to the subject). 

Proof. The fact that £& is independent of k for k > 1 follows from []TS] (pp 326 
ff). Choose k = 2. By 



2 a w 24^ w 6! 
Putting (|38a[) , (ggg ) and ( ggg) into (|§, we find 



OC' 



% (1) (0) + §^ (3) (0) - / Mt)9 i6) (t)dt . (46) 



o 



which leads to (^) by fl36|). The term 0(A 2 ) in fl4"T|) comes from making the 
change of variable t' = -j-t in the integral in ( |46|) and taking into account that ipk 
is 0(1). 
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What if we choose k = 1? By (|4*3"| ) and 
have the second term 



Si(0) = 0, but, in (|0[), we still 



1 Z" 00 

We use the recurrence ( IfTSfl , 13.2.13) 

= (2m + 1) {^ 2m + (-l)™- 1 ^} 

with m = 2, obtaining 



V>4 - £ 2 = 5^5 



(1) 



which we insert in 



getting 



24 L 5 y w 24 



9 {4 \t)dt 



(48) 

(49) 
(50) 

(51) 



Integration by parts in the first term on the r.h.s. of (|5l"D and use of (|38cj ) in the 
second term yield (using ip n (0) = 0) 



Si 



120 X + 



A further integration by parts using the recurrence 

brings fl5"2] ) to the form (|46|). We have thus proved 

5 2 /vr\3 



13.2.13) 



(52) 



(53) 



(54) 



for all k > 1 (the present argument is easily generalized). Thus, for parallel plates 
and Dirichlet b.c. the A-independent term in the asymptotic series ([§) is regulariza- 
tion independent and is the (3R, 0) sum of the divergent series (|36|). Neumann b.c. 
yield the same result. 

3 The Sphere 

The Casimir effect for b.c. on a sphere was first considered in the classic paper by 
Boyer [Tl[ and since it has been considered from diverse viewpoints: source theory 
j25| , multiple scattering | 26fl , dimensional dependence of the effect |27j as well as an 
improved mode summation method |1^, [HJ (see also [0). In |31[] it is shown how 



a natural subtraction method ensures convergence of the mode sum and in [32] RI 
for the ball has been proved (for a more detailed reference list see |H| and |22]| ). 

Here we will to reconsider the CE for a massless scalar field subjected to Dirichlet 
b.c. on a sphere in the light of the above developed theory. We will consider the 



12 



original sphere, of radius a, embedded in a concentric greater sphere of radius R> a. 
As it is well-known, for the sphere it is convenient to consider the inner and outer 
regions together in order to avoid logarithmic contributions for the CE. So, taking 
into account the (21 + l)-fold degeneracy of each eigenvalue, we have 



oo ^ oo 



(55) 



1=0 



n=l 



where u n i are the eigenfrequencies. The sum over n in (55) can be changed into an 
integral by using the Cauchy theorem |T^, |TJ], ^3]. Here we will follow [I4~| with an 
crucial difference: the cutoff functions used in [[14}] , while appropriate to treat the 
eletromagnetic field, do not render the integrals well-defined in the massless scalar 
field case, so that we will adoptf] 



C nl (A) 



-h(v/a+u nl ) 



(56) 



for the cutoff functions, which satisfies the normalization condition (TJ). Besides this, 
it is important to analyse the asymptotic behaviour of E vac based on more general 
cutoff functions. 

Then, we can rewrite (1551) as |T5|, [TJ] 



E v 



-. oo 



(57) 



1=0 



with (y = l + \) 



Qi = — e- Al/ / a Re e~ iv 

7T 



y exp{-iu-ye~ iip }— In /, {%vye~ ilp ) 
a ay 



dy , (58) 



where <p is an (small) angle which orientates the contour of integration with respect 
to the imaginary axis of z (see ||14|), and 



fi(iz) 



-IJz)KJz) 



(59) 



Now using the uniform asymptotic expansions for the Bessel functions l v and K v 



30| we can obtain an asymptotic expansion for Qi which is valid for large orders. 
Then, in general, we can rewritte E vac as (15 



E. 



asym 



1 - 



(60) 



1=0 



where ^asym stands for the expression obtained from (|57| ) - (|58|) by using the asymp- 



totic expansions for the Bessel functions |30[] , AQi = Qi — Qf 37 " 1 and n is such that 



1 



■i-More general cutoffs of the above type have been used by C. R. Hagen in a different context 
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for I > n the asymptotic expansion Q^ sym affords a good approximation for Qi (i.e., 
AQi ~ for I > n). 

Then, we obtain (after performing a rotation of the integration contour ye~ lip — ► 

y) 



E : 



asym 



2 a 2 

7T A 3 



2/ 



1 1 



7T 



A 



Re 



(1 +iy) 3 dy 
y d 



hit dy 



;i + %) dy 



C 2, - Re / y— 

ira 2 J dy 



a(£) dy 



(61) 



dy + 0(A) 



where ((s,a) = Y^o(^ + a )~ s * s the Hurwitz zeta function. From this expression 
must be clear why we have introduced the cutoff functions fl56|) rather than e~ AuJnl 
used in |TJ]]. Namely, in the absence of the term e~ kv l a in (M) the first integral 
in (31) would have a non-integrable singularity in the origin, but all integrals are 
well-defined if we adopt (|56|). 
From (H) we have (C(2, §) = 



35tt 5 



asym 



8 A 3 1024 A 65536 a 



+ 0(A) . 



(62) 



Now, it remains to calculate ^]™ =0 AQz in (|60|) . Notice that in this term the sum 
is finite and we do not have any divergence. Then, since ip > may be considered a 
small angle (sin<£> > and cosip > 0) we may integrate (|58f) by parts and perform 
a rotation of the integration contour {ye~ lip — > y) to obtain 



Qi = -~ [ dy\n[2yI„(y)K v (y)} + 0(A) 
* Jo 



(63) 



which is nothing but the Qi in ref. [15] (except for a sign). So we may take advantage 
of the numerical results in |15] for this expression. 

Analogously, we may obtain a expression for Q^ sym , appropriate for when there 
is no infinite summation, given by 

,2 



Qi 



asym 



V 
7T 



dy^— In 
dy 



y 



dy 



l 



dy a(t) 



0(t) - -a 2 (t) 



O(A) 



which after integration yields 

gasym = 



1 



2 + 128 



35 



32768 v 2 



+ 



(64) 



(65) 
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Then, we can take n = 4 in (pD|) as a good approximation (see flT5"| ) obtaining 



a 2 5 0.002819 , , 

Evac = -8A^- 1024A + ^^ + ---' (66) 

which yields a4 Phere ~ 0.002819 for the coefficient of the A-independent term in the 
asymptotic series (|3[) for E. Therefore the CE is 

af here 0.002819 , , 

^Casimir = " ^ , 67 

a a 

and by fl2"4j) we see that the Casimir force for massless scalar field with Dirichlet b.c. 
on a sphere is repulsive. This result was obtained with greater precision in |27| (also 
see 0,0). 

While the numerical result provided by ([67] ) is not new, the above calculation 
illustrates the fact that when we use more general cutoffs like fl5"6|) (which in the 
present case is mandatory) we are faced with an asymptotic series in A for E Yftc , see 
(|3|) and ( |66|) . Then, the method discussed above provides an unambiguous way to 
identify the CE. 



4 The Interior Problem for the Cube 



Consider now a cube K of side L, with Dirichlet b.c. (Neumann b.c. may be handled 
analogously). The normalized eigenfunctions and eigenvalues of (—A) 2 are 



3 

2V . /ni7ix 1 \ . (n 2 r nx 2 \ . fn 3 7rx 3 \ 

— sin sin sin 

L \ L J \ L J \ L J 



-A) 2-u r 



71 



n\u m 



n2n 3 {%) ] 



\n\ 



{nl+n 2 2 + nl)* 



where n« = 1, 2, 
We consider 



1,2,3). 



£ vac (A) = / d 3 xH(x,A) 



K 



By (H, 
By (|D, u 

^vac(A) = 



1 d 

£vac(A) = - — 



j\n\ and hence 



1 2tt 3 J V 

k. n 



2vr 2 16 dA 



.ngZ 3 



•a|n| 



3 e ' aVA +zJ2 e ~ al 

Hez 2 



»l 



(68) 



(69) 



(70) 



(71) 
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where 



71 A 

a = —A . 

Lj 



(72) 



The last sums in ( [TIP are due the fact that, because of (|68D , the planes n\ = 0, 
n2 = 0, and n3 = have to be excluded from the sum over Z 3 because they lead 
to eigenfunctions which are zero. For the same reason the axes n\ = n 2 = 0, 
n i — n 3 — 0, n 2 = n 3 = and the origin be excluded. Exclusion of the three planes 
(the term — 3 X^nez 2 e _a ' n ' in ([71])) corresponds to excluded each axis twice instead 
of only once. The third term compensates for this, while the last one excludes the 
origin. 

A method of calculation of the lattice sums in ([71] ) is through the Poisson sum- 
mation formula 

where are the Fourier coefficients of /: 



C n 



d x e 



-2-nim-x 



(74) 



See also ref. |0|. Applying ([751 ) to ([71]), we find 
3 



£vac(A) 



3 A -4 



2tt 2 



-LA 



3 L 3 A -4 _ A L 2 A -3 + A LA -2 

2vr 2 4tt 8tt 



mGZ 3 
m#0 



7rA 



+ 4ir 2 \m\ 2 



mez° 

m#0 



^- ) + J rr~ | /H 1 ~ 



8L \L 



7rA 



+ 4n 2 \m\ 2 



(75) 



— ) + \n 2 \lT) 1 ' 2 



3tt 
8L 



mez 



7T 



A 



L 



+ A7i 2 m 2 



AL \L 



m,eZ 



^— j +47T 2 m 2 



We now expand the sums X]m^o m dUD ^ n ^ ne following way: 

sA 2 



^— ) + 4ir 2 \rh\ 2 



(An 2 \m\ 2 y s (l 



4L 2 |m| 2 



+ 



(76) 



The unit term in (|76|) yields a contribution of type a±L 1 in (|3]), the remaining terms 
provide the rest of the asymptotic series in (Q) consisting of positive powers of A. 
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We thus find 

3 , 3 



04 



y |^|-4 + * y |^|-3 _ y m - 2 (77) 



32tt 2 ^ ' ' 64tt ^ ' ' 32tt 

mSZ 3 mGZ 2 
m^O 



The last sum above is nothing but 2£(2), where ( stands for the Riemann zeta 
function, and the second one may be rewritten as the product of two independent 
sums by means of X^ez 2 l^l -5 = (f) P (§) (see, e.g., [20, |23|), where /3(s) = 



Yl'jLo (27T1T 7 ' Then, using the result of Lukosz [^4| for the first sum in (f7?D we 
obtain 

a 4 = -0.0157322... , (78) 

which is in accordance with the result obtained numerically in ref. [^D] (in fact we 
have obtained 04 to a higher accuracy than shown). In addition, from (|23b ), the 
inner pressure is 

pJTW = . (79) 

By (|75|) and (|79|) we see that the force due to the interior is attractive. The 
repulsive character of the sphere (section ^) suggests, however, that the same is true 
for the cube. This fact alone shows that this sign, if true, must be entirely due to 
the exterior force, a subtle problem to which we now turn. 



5 The External Problem for the Cube 

As remarked above, it is of great interest to consider also the outer problem for the 
cube. We will consider the cube Kl of side L concentric with a cube Km, of side 
M, from which Kl is a subset (M > L and M eventually goes to infinity at the end 
of calculation) and impose Dirichlet b.c. on Kl as well as Km (see section 1). Un- 
fortunately, the solution of the external Casimir problem for the cube with Dirichlet 
b.c. cannot be constructed out of the functions of the form Q58|), because the con- 
tinuity conditions on several planes cannot be satisfied simultaneously. However, 
the form of solutions which are naturally adapted to the internal geometry of 
the cube, suggest splitting the region Km\Kl into 26 subregions bounded by the 
planes containing the faces of the cube. We may require the u n (x) to vanish on the 
boundaries of these subregions, including the original requirement of vanishing on 
the faces of the internal and external cubes. If we do so, the resulting problem is 
explicitly solvable in terms of the set fl6"B|) . Of course, this Ansatz introduces addi- 
tional stresses in the region Km\Kl- We shall comment on these restrictions at the 
end of this section. 

Then we have that the 26 subregions which compose Km\Kl are of three topo- 
logically distinct kinds (with both cubes centered in the origin): 1) a rectangular 
box with two sides L and one M ~ - with multiplicity 6; 2) a rectangular box with 
two sides M ~ L and one L (with multiplicity 12) - the contribution of the edges; 3) 
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a cube of sides M ~ L (with multiplicity 8) - the contribution of the corners. The 
Casimir energy of each of these regions can be obtained along the same lines of the 
calculation above outlined for the inner cube (see JTS|). Then we obtain that the 
regions of type 3) do not contribute, i.e., there are no contributions of the corners, 
either to the energy or the pressure, in the limit M — > oo. The total contribution 
of the regions of type 1) is 



E\(L, M) 



3L 2 (M-L) 
32^ 



E 



mSZ° 



3L(M - L) j2 



m\L 2 + m|L 2 + m| < ^-r — 



1 



32tt 



•»sz 2 m? 



L 2 + 



TO 



2 (M-L) 2 



+ 



32ttL 



E 



2 4 
77 / 2 



(80) 



2 [mf + m\\ 



16 V L + M 



and for the regions of type 2) the total contribution is 

3L(M-L) 2 



E 2 (L : M) 



+ 



+ 



327T 2 

3L(M-L) 
3 



mez 3 



1 



?L 2 + m 



2 (M-L) 2 



mGZ 2 m? 



L 2 + m\ 



2 (M-L) 2 



(81) 



8tt(M - L) 



rnE 



2 L m^ + m|] 2 



7T / 1 
3 - - I — + 



4 



8 \L M-L 



Then, by ( |23c|) and taking into account that dV onte r = 6dVi + YLdVi + 8dVz, where 
Vi is the volume of type i) region, and that in the thermodynamic limit dV$ — > 
(dV 3 /V 3 oc 1/M), we obtain 

p° uter (L) = - lim — ^ -A M) + £ 2 (L, M)l . (82) 

yc V ; m^oo 3M(M - 2L) 1 V ; 

It may also be verified explicitly that regions of type 3) do not contribute to the 



energy, or to the pressure, in the thermodynamic limit. From (|82|) we see that type 
1) regions do not contribute to the outer pressure, while the edges contribution is 
given only by the first term in (31 

P°J tcr (L) 



( 



lim \ ^~ 



L 



32tt 2 



V 



m 3 ^0 



m\L 2 + m| 



2 (M-L) 2 2 (M-L) 2 



r 

/ 



(83) 
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It follows from 



that (see later): 
/ 



32tt 2 



£ 

m#0 



/ 



dL L 



7T~ 



480 



-L" 



(84) 



There is no contribution in the thermodynamic limit to ([83]) from the sum over 
m2 7^ or m 3 ^ 0, or both. Indeed, the term m\ = in (|83| ) does not contribute as 
M — > oo, as one sees easily, and 



mj^O L 



2 (M—L) 2 



+m 



2 (M-L) s 



n 2 



< 



f 2 

\L 7 (M-L) 



£- 



1 



(85) 



If only r?72 7^ (or 7773 7^ 0) 



£ 



|mf L 2 + 777,2 



< 



2 (M-L)2 



L 3 (M-L) 



E 



777,eZ^ 
7^0,7712 7^" 



7?7| + 777r 



(86) 



for L < 1, M > L. It may be checked that these bounds suffice to show that the 
contributions of the above sums to the limit in ( |8"B| ) is zero. The only surviving term 
in the sum in (T~ 



_J is thus 7?72 = 7773 = 0, which leads to 
It is most important to note that the edges' contribution to the Casimir pressure 
is greater than the inner pressure in absolute value: by 



3) and (pal) 



p c (L) = (-0.005244 + ^L" 4 



0.015317 L~ 



(87) 



The Casimir pressure is thus repulsive, and the net result is that edge effects deter- 
mine the sign of the force. We shall return to this point in the conclusion. 

We now comment on our Ansatz. We have replaced the original Hilbert space 
of L 2 -functions with Dirichlet b.c. in the inner and outer boundaries with a direct 
sum of 26 spaces, upon introduction of additional Dirichlet b.c. on planes which are 
extensions of the cube's faces to the region Km\Kl- In this region all extra stresses 
are parallel to the cube's faces, and for this reason the Casimir pressure is insensitive 
to their inclusion. We have verified this assertion in appendix A by introducing and 
extra Dirichlet plane orthogonal to a system of parallel plates in the inner region. 
The proof generalizes to an arbitrary finite number of such planes provided they are 
placed symmetrically to some plane orthogonal to the z-axis, thus not introducing 
an extraneous length in the original problem. Since the parallel plates are a soluble 
"limiting case" of the cube, we (strongly) believe that our Ansatz provides the exact 
solution for the cube. 
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6 Conclusion and Open Problems 



In this paper we have introduced a mathematically precise framework for the Casimir 
pressure, by associating it to the cutoff-independent part of the Ramanujan sum of 
the (divergent) series for the Casimir energy. Our ideas have precursors in |l0] and 
|§. We illustrated the framework by parallel plates, the sphere and the interior 
problem of the cube. 

In section [5] we introduced an Ansatz to calculate the exterior Casimir pressure 
for the cube. We discussed why we (strongly) believe that it is the exact solution 
for the cube. If our conjecture is right, the calculation of section |5| also provides an 
explanation for the sign of the force: it is due to a competition betweeen the inner 
and outer pressures, in which the latter is positive and larger than the former in 
absolute value, because, as remarked in section the thermodynamic limit selects a 
set of modes different from the inner ones|, with a large positive contribution from 
the edges. The edges reflect the passage from the infinitely extended parallel plates 
to a compact region, i.e., by folding. If this folding were smooth, i.e., for any smooth 
approximation to the cube, it would be accompanied by nonzero curvature. At the 
other extreme - uniform nonzero curvature - we have the sphere. Here, however, 
curvature effects appear less directly, reflecting themselves in the appearance of the 
Neumann functions in the external problem. It is an interesting open problem to 
understand more clearly the role of curvature (of various kinds, e.g. Riemannian, the 
mean and Gaussian curvatures) in the Casimir effect for general compact manifolds 
with boundary (see also fT0|l ). 
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Appendix A 

In this appendix we consider the problem of the parallel plates with an orthogonal 
(Dirichlet) plane introduced in the inner region and prove that the Casimir pressure 
on the plates is the same that for the problem without the Dirichlet plane (see 
section g). 

Consider the two parallel plates placed at z = and z = d, and the inner Dirichlet 
plane at y = 0, < z < d. The outer problem is the same which we worked out in 

§See also [[33] for a discussion (different of ours) of the different roles of the inner and outer 
modes of the Casimir problem 
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section |2|, and it does not contribute. The inner Casimir problem is now split into 
two regions: fq nner = {x G Ix [0, oo) x [0, d]} and K' 2 nner = {x G Ex(-oo, 0] x [0, d]}. 
In order to calculate the Casimir energy of these regions it is necessary to consider a 
finite region in the plane x-y with area A = L\L 2 (the whole area of the plates), and 
take the limit L%,L 2 — *■ 00 at the end. Then, both K™ ner and K^ ner are given by 
boxes of sides Li (along x-axis), L 2 /2 (along y-axis) and d, so that we can proceed 
in the same way that in section f|, obtaining for the A-independent term of the series 
d) for the region K[ nnei : 



E™ er (L u L 2 ,d) 



LiL 2 d 
64vr 2 



E 



1 



+ 



L\L 2 

1287T 



L 2 d 

1287T 



m\L\ + m 2 .^ 2 - + mid 2 



Lid \ — ■> 



y l 

t t \m\L\ + m?4 1 § 6471 ^ K L i + 



E ! — 



7T 



1 2 1 
96 \ ^1 + + d 



, (A.l) 



and, obviously, E 2 nnev (Li, L 2 , d) has the same form. 

Now we can calculate the inner pressure by means of ( |23b|) , where V\ T 
be taken as the whole interior volume: Vinncr = L\L 2 d. Thus, we have 



must 



Pc{d) 



inner 

Pc 



(d) 



lim 

L\ ,L,2— >oo 



lim 

2_d_ 

LiL 2 dd 



dET er (Li,L 2 ,d) 
8V 

u * inner 

E^ nex (Lx, L 2 , d) 



8V 

v * inner 



(A.2) 



where we have taken into account that E™ ner = E 2 nncr . 

From ( |A.1J ) and (|A.2|) we see that the only term which contributes to the pressure 
at the thermodynamic limit is the first term at the right side of ( |A.1[) ; in fact, only 
the term ni\ = m 2 = contributes (that the terms with mi 7^ or m 2 7^ 0, or both, 
do not contribute can be proved just in the same way that in section || - see eqs. 
BT)). Then, we obtain 



Pc(d) 



I67H ad 



71' 



480 



-d~ 



(A.3) 



which is the same result obtained in the case without the inner Dirichlet plane! 

The above result shows that the introduction of additional stresses parallel to 
the physical plates do not modify the pressure on these plates. Finally, it is easy 
to see that the above proof generalizes in a trivial way to the case in which the 
Dirichlet plane is at the outer region, and also to the case in which we have a finite 
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number of Dirichlet planes in the inner region provided these planes are disposed 
symmetrically with respect to some plane orthogonal to the z axis. 



Note added. The following related references have been brought to our atten- 
tion: M. Bordag, E. Elizalde, K. Kirsten, S. Leseduarte, Phys. Rev. D 56, 4896 
(1997) for spherical geometries; and N. F. Svaiter, B. F. Svaiter, J. Phys. A 25, 979 
(1992) and Phys. Rev. D 47, 4581 (1993) for the use of auxiliary cavities in order 
to treat surface divergencies. 
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